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ABSTRACT 

We present the results of numerical experiments, in which we study how the 
asphericities induced by the growth of the standing accretion shock instability 
(SASI) produce the gravitational waveforms in the postbounce phase of core- 
collapse supernovae. To obtain the neutrino-driven explosions, we parameterize 
the neutrino fluxes emitted from the central protoneutron star and approximate 
the neutrino transfer by a light-bulb scheme. We find that the waveforms due to 
the anisotropic neutrino emissions show the monotonic increase with time, whose 
amplitudes are up to two order-of-magnitudes larger than the ones from the con- 
vective matter motions outside the protoneutron stars. We point out that the 
amplitudes begin to become larger when the growth of the SASI enters the non- 
linear phase, in which the deformation of the shocks and the neutrino anisotropy 
become large. From the spectrum analysis of the waveforms, we find that the am- 
plitudes from the neutrinos are dominant over the ones from the matter motions 
at the frequency below ~ 100 Hz, which are suggested to be within the detection 
limits of the detectors in the next generation such as LCGT and the advanced 
LIGO for a supernova at 10 kpc. As a contribution to the gravitational wave 
background, we show that the amplitudes from this source could be larger at the 
frequency above ~ 1 Hz than the primordial gravitational wave backgrounds, but 
unfortunately, invisible to the proposed space-based detectors. 
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Subject headings: supernovae: collapse — neutrinos — hydrodynamics — insta- 
bility 

1. Introduction 

The gravitational astronomy is now becoming a reality. In fact, the ground-based laser 
interferometers such as TAMA300 (Ando & the TAMA collaboration 2002; Ando et al. 
2005) and the first LIGO (Thorne 1995; Abbott et al. 2005) are beginning to take data at 
sensitivities where astrophysical events are predicted. For the detectors including GEO600 
and VIRGO, core-collapse supernovae especially in our Galaxy, have been supposed to be 
the most plausible sources of gravitational waves (see, for example, New (2003); Kotake et al. 
(2006) for review). Since the gravitational wave (plus neutrinos) is the only tool which gives 
us the information in the innermost part of evolved massive stars, the detection is important 
not only for the direct confirmation of gravitational waves but also for the understanding of 
supernova physics itself. 

So far, most of the theoretical predictions of gravitational waves from supernovae have 
focused on the bounce signal in the context of rotational (Monchmeyer et al. 1991; Yamada 
& Sato 1995; Zwerger & Mueller 1997; Dimmelmeier et al. 2002; Fryer et al. 2002; Kotake 
et al. 2003; Shibata & Sekiguchi 2004; Ott et al. 2004) and magnet orot at ional (Kotake et 
al. 2004; Obergaulinger et al. 2006) core collapse. In most of the previous studies, the iron 
core prior to core-collapse was assumed to rotate much more rapidly than predicted by the 
recent stellar evolution calculations (Heger et al. 2005). Recently, the initial rotation periods 
were estimated to be larger than ~ 100 sec for the observed rotation periods of the radio 
pulsars (Ott et al. 2006). In such a slowly rotating case, the bounce signal becomes too 
small to be detected even by the laser interferometers in the next generation for a galactic 
supernova, owing to the suppression of the rotation-induced deformation at core-bounce (see, 
e.g., Kotake et al. (2004)). 

Besides the rapid rotation of the cores, two other ingredients have been considered 
to be important in the much later phases after core bounce, namely convective motions 
and anisotropic neutrino emissions. Both of them contribute to the non-spherical parts in 
the energy momentum tensor of the Einstein equations, thus being the potential sources 
of the gravitational wave (see, Kotake et al. (2006) for a review). One of the possibility 
as the origin of the asphericities may be large scale density inhomogeneities formed in the 
central core prior to collapse (e.g., Bazan & Arnett (1998); Meakin & Arnett (2006)). Fryer 
(2004) performed three dimensional SPH simulations and pointed out that the neutrino- 
originated gravitational waves, which dominate over the one from the convections, are within 
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the detection limits for the advanced LIGO for the galactic supernova (see also, Burrows & 
Hayes (1996); Fryer (2004); Fryer et al. (2004); Muller k Janka (1997)). Another possibility 
to induce anisotropy is the (moderate) rotation of the core. Muller et al. (2004) calculated the 
gravitational waves based on the two-dimensional (2D) Boltzmann transport simulations of 
slowly rotating core (Buras et al. 2003) and found that the neutrino-originated gravitational 
waves exceed the bounce signal large enough to be detectable by the advanced LIGO with 
good signal-to-noise ratio for the galactic supernova (see, also Kotake et al. (2006) for the 
properties of neutrino-originated gravitational waves in the rapidly rotating case). More 
recently, the new ingredient of the gravitational- wave emissions is reported (Ott et al. 2006), 
namely the g-mode excitations of the protoneutron stars, which was observed in the 2D 
approximate Boltzmann transport simulations at much later postbounce phase (~ 600 ms) 
(Burrows et al. 2006). 

There is an another ingredient for producing large asphericity, to which much attention 
has been paid recently in the context of the studies about the explosion mechanisms, that 
is the so-called standing accretion shock instability (often called "SASI"). In the numerical 
simulations by Blondin et al. (2003); Scheck et al. (2004); Blondin & Mezzacappa (2006); 
Ohnishi et al. (2006a,b), it was found that the standing shock wave is shown to be unstable to 
non-radial perturbations, and that the perturbations grow up to the non linear regime with 
clear low-mode (£ = 1, 2) dominance, leading to the global deformation of the shock wave 
later. Here £ stands for the azimuthal index of the Legendre polynomials. The importance of 
SASI is also stressed by the recent studies, demonstrating that such an explosion is favorable 
to reproduce the observed synthesized elements of SN1987A (Kifonidis et al. 2006) and also 
to explain the origin of the natal kicks of young pulsars (Scheck et al. 2004). These situations 
motivate us to study how the gravitational waveforms are originated from the asphericities 
by SASI. 

In this paper, we present the results of numerical experiments, in which we study how 
the asphericities induced by the growth of SASI produce the gravitational waveforms. To 
obtain the neutrino-driven explosions, we parameterize the neutrino fluxes emitted from the 
central protoneutron star and approximate the neutrino transfer by the light-bulb scheme. 
Based on the long-term two dimensional hydrodynamic results, we calculate the gravitational 
waveforms. By doing the spectrum analysis, we study the detectability of such signals from 
a nearby core-collapse supernova. It is noted that much attention has been paid recently 
to the core-collapse supernovae as one of the promising sources of the cosmological gravita- 
tional wave backgrounds (see, e.g., Buonanno et al. (2005) and references therein). Thus we 
calculate the SASI-induced gravitational wave backgrounds and discuss the detectability by 
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the currently proposed space-based detectors such as LISA 1 , BBO 2 , and DECIGO (Seto et 
al. 2001). 

The plan of this paper is as follows: In Section 2, we outline the initial models, the 
numerical methods, and shortly summarize the methods for calculating the waveforms. We 
show the main numerical results in Section 3. We summarize and discuss our results in 
Section 4. 



2. Numerical Methods and Models 

2.1. Initial Models and Hydrodynamics 

The numerical methods employed in this paper are essentially the same as those used 
in our previous paper (Ohnishi et al. 2006a). The basic evolution equations, describing the 
compressible accretion flows of matter attracted by the protoneutron star and irradiated by 
neutrinos emitted from the neutrino sphere, are written as follows, 

f t+P V-v = 0, (1) 



^dt ~~ ' V + Q E + Qincl — Qvu — Qbrcms ~ Qplasmon, (3) 

f = Q„ (4) 

4>=-^, (5) 
r 

where p, v, e, P, Y e , $ are density, velocity, internal energy, pressure, electron fraction, and 
gravitational potential, respectively. We denote the Lagrangian derivative as d/dt and r is 
the radius. Mj n is the mass of the central object. The self-gravity of matter in the accretion 
flow is ignored (see Yamasaki & Yamada (2006) for the effect). The parameters of Qe and 
Qn are related to the standard heating and cooling via neutrino absorptions and emissions 
by free nucleons (see also Ohnishi et al. (2006a)). Q mc \ is the minor additional heating by the 
inelastic neutrino- helium interactions estimated by Haxton (1988) as considered in Ohnishi 



^ee http://lisa.jpl.nasa.gov/ 

2 see http:/ /universe. nasa.gov/program/bbo. html 
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et al. (2006b). In addition, we newly take into account the neutrino cooling by neutrino pair 
annihilation to e~e + pairs (: Q u p), nucleons-nucleons bremsstrahlung (: Qbrcms) (Hannestad 
& Raffelt 1998; Burrows et al. 2000), and plasmon decay into u e , u e (: Qpiasmon) (Schinder et 
al. 1987; Ruffert et al. 1996). 

The numerical code employed in this paper is based on the modified version of the ZEUS- 
2D (Stone & Norman 1992) for the applications to the supernova studies (Kotake et al. 2003), 
in which the tabulated realistic equation of state (EOS) based on the relativistic mean field 
theory (Shen et al. 1998) was implemented. Furthermore, we have added the equation for 
electron fraction (Eq. (4)) and included the neutrino coolings/heatings parametrically as the 
source term of the energy equations (Eq. (3)), both of which are solved in the operator- 
splitting fashion. In the simulations, spherical coordinates are used without imposing the 
equatorial symmetry. The computation domain covers the whole meridian section with 60 
angular mesh points, except for a model in which we have adopted 120 angular grid points. 
Since the latter model did not produce any significant difference from other models, we will 
report in the following the results obtained from the models with 60 angular mesh points. 
We use 300 radial mesh points to cover r m < r < r ont = 2000 km, where r m is the inner 
boundary and chosen to be roughly the radius of neutrino sphere. 

The initial conditions are provided in the same manner of Ohnishi et al. (2006a) as the 
steady state solution of Yamasaki & Yamada (2005). In constructing the initial conditions, 
we assume a fixed density p- m = 10 11 g cm -3 at the inner boundary. And the initial mass 
accretion rates and the initial mass of the central object are set to be M = 1 M s" 1 and 
Mi n = 1.4 Mq, respectively. To induce the non-spherical instability, we have added £ = 1 
velocity perturbations to the initial state mentioned above. At the outer boundary, we 
adopt the fixed boundary condition consistent with the initial condition. On the other hand, 
the absorbing boundary is used at the inner boundary. The temperatures of electron-type 
neutrinos are also constant and set to be T Vc = 4 MeV and T Pc = 5 MeV, which are the 
typical values in the post-bounce phase. The temperature of mu and tau neutrinos is chosen 
to be T v = 10 MeV. In the standard model, the luminosity of electron-type neutrino L Uc 
and anti-neutrino L Pc are set to be 6.5 x 10 52 ergs s -1 . In addition, we examined two cases 
of lower luminosities of L Uc = L Pc = 6.0,5.5 x 10 52 ergs s _1 . In all the computed models, 
the luminosity of mu and tau neutrinos is set to be half value of electron-type neutrinos, 
keeping consistency with the results obtained by the previous detailed numerical studies 
(e.g., Liebendorfer et al. (2001)). 
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2.2. Computations of gravitational wave signals 



We follow the methods based on Epstein (1978); Miiller & Janka (1997); Burrows & 
Hayes (1996); Kotake et al. (2006) in order to compute the gravitational waveform from 
anisotropic neutrino emissions. We will summarize it in the following for convenience. Since 
we assume axisymmetry, the transverse-traceless gravitational field h TT is shown to have one 
nonvanishing component for the observer in the equatorial plane: 



h 



TT 



AG 



dt / dte' $(0') 



dl v (0',t') 



dn' 



(6) 



where G is the gravitational constant, c is the speed of light, R is the distance of the source 
to the observer , dl v /dVL represents the direction-dependent neutrino luminosity emitted per 
unit of solid angle into direction of f2, and denotes the quantity, which depends on the 
angle measured from the symmetry axis (#') (see Figure 1), 



$(0' 



smi 



-TT 



2n 



d(f> 



1 + sin 9 cos i 



, Q 1 + tan 2 & sin 2 0' 
7rsin0'(-l +2|cos0'|). 



(7) 
(8) 



The integration with respect to the azimuthal angle (0'), albeit fairly straightforward, is 
pretty helpful to the later discussions, which is not published to our knowledge. In addition, 
from Eq. (8), one can readily see that no gravitational waves are emitted if the neutrino 
radiation is isotropic. We estimate dl u /dVL as follows, 



dl v 

dn 



(0 



dr r 2 \Q C \ 



calculated at each $' 



(9) 



where Qq is the sum of contributions from the neutrino coolings, of each species (u e , is e , 
and ^u, r ) computed from the source term of Eq. (3). In the above estimation, neutrinos 



are assumed to be emitted radially in each angular bin. As a guide to see the anisotropy of 
the neutrino emissions, we calculate the anisotropy parameter according to Miiller & Janka 
(1997) as follows, 



wo 



It should be noted that the gravitational waves from neutrinos have a different feature from 
the ones from matter motions, in the sense that the former has memory effect, which means 
that the gravitational amplitude jumps from zero to a nonvanishing value and it keeps the 
non-vanishing value even after the energy source of gravitational waves disappeared (see 
Braginskii & Thorne (1987) for details, and Segalis & Ori (2001); Hiramatsu et al. (2005) 
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Fig. I. — Angular dependence of <3> in Eq. (8). Note that the angle is measured from the 
symmetry axis. 



for the examples of the astrophysical emitters of such gravitational waves and references 
therein). In Eq. (6), this nature can be directly seen as the time-integral. 

As for the gravitational waves, hit), of the quadrupole radiation of mass motions, we 
employ the standard strain formula, 

1 /1 \ 1/2 A E2 

where the form of A^q can be found in references (see equation (12) in Monchmeyer et al. 
(1991)). In using the formula, the azimuthal gradient of the gravitational potential is set to 
be zero, since we neglect the self-gravity of the computational domain and treat the gravity 
as in Eq. (5). In the following computations, we assume that the observer is located in the 
equatorial plane since the most of gravitational wave is radiated in the plane (a = tt/2 in Eq. 
(11)), and that the source is assumed to be located at our galactic center (R = 10 kpc). In 
order to assess the detectability of the gravitational waves obtained in this study, we employ 
the characteristic gravitational wave strain, 



h m - 1 / 2 G dE ow(f) ri9 v 

M/) -i?V^^/^' (12) 

for a given frequency / (Flanagan & Hughes 1998). Here dEcw(f)/df is the energy spectra 
of gravitational waves defined as follows, 

dE ow (f) _ . 



df G 16tt 
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We employ the Fast Fourier Transformation technique in order to perform the Fourier trans- 
formation of A^{t) to if 2 (/)- 

We discuss the contribution of the gravitational waves mentioned above to the cosmo- 
logical background gravitational radiation. For simplicity, we assume that all core-collapse 
supernovae have identical emission characteristics. According to Phinney (2001); Buonanno 
et al. (2005), the density parameter, which represents the sum of energy densities radiated 
by a large number of independent core-collapse supernovae at each redshift, can be written 
as, 



1 f 00 
flgw(/) = — ^ / dz 

Pc<~- Jo 



Rsn(z) 



dt 



dz 



l + z 

where p c = 3Hq/(8tcG) is the cosmic critical density, Rsn(z) is the event rate of core-collapse 
supernovae per comoving volume, and f z = (1 + z)f. The cosmological model enters with 
\dt/dz\ = [(1 + z)H(z)]' 1 and, for a flat geometry, 

H(z) = H [n M (l + zf + n A ] 1/2 . (15) 

We will use the parameters Qm = 0.3, VIa = 0.7, and H = ho 100 km s~ x Mpc -1 with 
h = 0.72. As for the core-collapse supernova rate, we employ the parameterization by 
Buonanno et al. (2005) as follows, 

R SN (z) - R SN x | 2/3 _ Q(i + z)a {oil < z < 2Q ■ (16) 

We take (3 = 2.7 and the present-day rate i?g N = 2 x 10~ 4 Mpc~ 3 yr _1 , which are 
consistent with the Super-K limits on the diffuse neutrino background (Malek et al. 2003). 
Since the parameter < a < 2 is much less constrained than (3 and i?g N (see, for example, 
Ando & Sato (2004) and references therein), we examine the two cases of a = 0, 2 in the 
later discussion. 



3. Results 

Only in the case of L Ve = 6.5 x 10 52 erg s -1 , we can observe the continuous increase of 
the average shock radius, reaching the outer boundary of the computational domain (2000 
km in radius) in ~ 500 ms with the explosion energy of 4.0 x 10 50 ergs. In other models, 
we terminate the simulations at about 1 sec, not seeing the increase of the shock radius. 
Bearing in mind the evidences that the SASI-induced explosions are favorable for explaining 
the observed quantities of supernovae mentioned earlier, we take the exploding model of 
L u = 6.5 x 10 52 erg s' 1 as a reference in the following. For later convenience, the values of 
several important quantities are summarized in Table 1. 
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Fig. 2. — Gravitational waveforms from the sum of the anisotropic neutrino emissions and 
the matter motions (left panel) and only from the matter motions (right panel). The time 
is measured from the epoch when the neutrino luminosity is injected from the surface of the 
neutrino sphere. Note that the supernova is assumed to be located at the distance of 10 kpc. 

3.1. Properties of Waveform 

In the left panel of Figure 2, the total amplitudes including the contributions both from 
the anisotropic neutrino emissions and the matter motions, are shown. Comparing the right 
panel, in which the amplitudes only from the matter motions are shown, we can see the 
neutrino-originated gravitational waves are up to two order-of-magnitudes larger than those 
from the matter motions during the SASI operation (at t > 200 ms as explained shortly). 
This reflects the small mass in the regions outside the neutrino sphere in the postbounce 



Table 1. Model Summary 



L v (10 52 crg/s) 


At (ms) 




/eq (Hz) 






6.5 


540 


6.1 x 10- 10 


51.0 


5.3 x 10" 


22 


6.0 


1000 


4.6 x 10~ 10 


59.0 


7.0 x 10~ 


22 


5.5 


1000 


2.3 x 10 -10 


63.2 


5.1 x 10- 


22 



Note. — L v denotes the input luminosity. At represents the simulation 
time. Egw is the total radiated energy in the form of gravitational waves in 
unit of Mqc 2 . / eq and h C q represents the frequency, and amplitude, under 
which the gravitational waves from anisotropic neutrino emissions dominate 
the ones from the matter motions, respectively. Note that the supernova is 
assumed to be located at the distance of 10 kpc. 
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phase (at most O.1M throughout the simulation time). In the following, we pay attention 
to the neutrino-originated gravitational wave. 

Looking at the left panel of Figure 2, the waveforms show the monotonic increase with 
time regardless of the input neutrino luminosities. To understand this trend, we look into the 
hydrodynamical behaviors induced by the SASI, taking the case of L Ue = 6.5 x 10 52 erg s _1 . 

In Figure 4, snapshots showing the hydrodynamical features with the resulting gravita- 
tional wave amplitudes (inserted figures) are shown. Up to ~ 100 ms after the onset of the 
simulation, no significant changes in the amplitudes and no deviations of the dynamics from 
spherical symmetry are found. It is noted that no gravitational waves are emitted when the 
motion of the regions outside is spherical, since the neutrino luminosity from the center is 
taken to be isotropic. 

After ~ 100 ms, the regions outside the neutrino sphere begin to oscillate with increasing 
average radius and wave amplitudes. The reason why the sign of the growing amplitudes is 
positive is as follows. From the angular dependence of <3> in Figure 1 and Eq. (8), it can 
be seen that the neutrino emissions from the regions with 9' smaller than n/3 contribute 
to the positive amplitude, while the negative sign comes from the regions from 7r/3 to tt/2. 
It is noted that this feature is north-south symmetric ( Eq. (8)). Due to the dominance of 
I — 1, 2 modes of the deformed shock waves during the SASI operation mentioned below, 
the neutrino emissions become stronger in the regions close to the symmetry axis. By these 
two factors, the amplitudes are found to increase monotonically with time. This property is 
common to the lower luminosity models as seen from the left panel of Figure 2. As a side- 
remark, the contribution of each neutrino species to the waveform is presented in Figure 3, 
which also shows the dominant contribution from v e . 

Next we discuss in more detail how the features of the waveforms are related to the 
growth of the SASI. For the purpose, we decompose the fluctuations of the shock surface 
into the spherical harmonic components; 



£=0 

Since the system is axisymmetric, only m = harmonics, that is Legendre polynomials, show 
up. The coefficients, ag, can be calculated by the orthogonality of the Legendre polynomials; 




(17) 




(18) 



The position of the shock surface, R s (9), is estimated from the iso-entropic surface of s = 5. 
The left panel of Figure 5 shows the evolution of the amplitude of each mode (£ = 1,2) 
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Fig. 3. — Contribution of each neutrino species to the waveform for the case of L Ve = 
6.5 x 10 52 erg s _1 . Here u x represents /i, r neutrinos. It is noted that the qualitative features 
are common to the other cases. 

normalized by the average shock radius of ao- The dominance of the fundamental mode 
{£ = 1) initially is simply due to the initial velocity perturbation assumed in this simulation. 
From ~ 40 ms, £ = 2 begins to develop rapidly and at ~ 200 ms, the amplitude becomes of 
the same order as that of the fundamental mode, which has already been saturated by this 
time. This marks the beginning of the nonlinear phase. 

This transition from linear to nonlinear phase corresponds to the time of the rapid 
increase of the average shock radius as seen from the top right (~ 200 ms) to middle left 
panel (~ 270 ms) in Figure 4. Simultaneously, this makes the deformation of the shock 
more elongated along the symmetry axis, leading to the slightly steep rise of the neutrino- 
originated gravitational wave afterwards for the reason mentioned above (see the left panel 
of Figure 2). Here it should be noted that the dominance of £ = 2 mode induces the rapid 
increase of the gravitational waves from the matter motions, albeit its amplitudes being 
very small, and is clearly understood from the quadrupole nature of gravitational radiations 
(compare the left and right panels of Figure 5). 

The anisotropy parameter a(t) in Eq.(10) also helps us to see the relation between the 
anisotropy of the neutrino radiation fields and the properties of the waveforms. The time 
evolution of a is presented in Figure 6. In case of L Ve = 6.5 x 10 52 erg s -1 , a takes larger values 
after the saturation, leading to more greater wave amplitudes. In the right panel of Figure 
6, the ratios of the wave amplitudes of L Ve = 6.5 x 10 52 erg s _1 and L Ue = 6.0 x 10 52 erg s _1 
to L Ue = 5.5 x 10 52 erg s _1 , are shown. Only near after the nonlinear phase sets in (~ 200 
ms), it is seen that the SASI-induced anisotropy of the neutrino emissions determines the 
waveforms, whose amplitudes monotonically grow with time. As a reference, we draw the 
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Fig. 4. — Entropy- (the left half panel of each panel) and net neutrino emissivity- (the right 
half) distributions in the meridian section for the case of L Ve = 6.5 x 10 52 erg s" 1 . The insert 
of each panel shows the wave amplitudes, in which the green point indicates the time of the 
snapshot. 
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lines indicating the ratio of the input neutrino luminosities. It can be also seen that the 
increase of the amplitudes is greater than the one estimated only by the difference of the 
input luminosities after the onset of the nonlinear phase. 

Finally, we have to discuss whether the above discussions are subject to change when 
the numerical resolution becomes more better. In Figure 7, we compare the waveforms when 
the angular resolution is doubled in comparison with the model computed so far. As easily 
guessed, the difference appears after the SASI saturates (~ 200 ms) and the non-linear phase 
sets in. In the lower resolution calculation, the anisotropy of the shock propagation tends 
to become more larger than the high resolution case, which results in the larger neutrino 
anisotropy, and thus leading to the larger wave amplitudes. The difference is not greater 
than ~ 20%, and thus the qualitative features discussed so far are found to be unchanged. 




Time [ms] Time [ms] 

Fig. 5. — Time evolution of the normalized amplitudes of I = 1,2 modes (left) and the 
waveform from the matter motions (right) for the case of L Ve = 6.5 x 10 52 erg s -1 . 



3.2. Spectrum Analysis 

Now we move on to discuss the features of the waveforms by the spectrum analysis. 
From Figure 8, one can see the dominance of the neutrino-originated gravitational waves 
at lower frequency of ~ 100 Hz. As seen in the left panel of Figure 2, the waveform from 
the neutrino shows a long-time variability in comparison with the rapidly varying (O(ms)) 
waveforms from matter motions, which is due to the local hydrodynamical instabilities. It 
is the memory effect of the neutrino-originated gravitational waves mentioned earlier that 
absorbs the rapid time variations of the neutrino anisotropy. 

To see clearly the dominance of the neutrino-originated gravitational wave than the 
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Time [ms] Time [ms] 

Fig. 6. — Time evolution of the neutrino anisotropy parameter (see Eq. (10)) (left) and 
the ratios of the wave amplitudes normalized by the case of the lowest input luminosity of 
L Ve = 5.5 x 10 52 erg s _1 (denoted as L 55 ) (right). The horizontal lines in the right panel just 
represent the constant ratio of the input luminosity. 

matter-originated one at the lower frequency, we define / eq , the frequency below which the 
dominance occurs, and the corresponding gravitational wave amplitude, h eq (see Table 1). 
From the table, it is found that the typical frequency is smaller for the higher luminosity 
case. In Figure 9, the gravitational wave spectra are plotted with the sensitivity curves of 
the laser interferometers. It can be seen that the detection of the gravitational wave at the 
low frequency range becomes more promising thanks to the contributions from the neutrinos. 
It can be seen that the gravitational waves from the neutrinos, which are dominated below 
~ 100 Hz, seem marginally within the detection limits of the currently running detector of 
the first LIGO and the detection seems more feasible for the detectors in the next generation 
such as LCGT and the advanced LIGO if a supernova occurs in our galactic center. 

Finally, we are in a position to discuss the contribution of the gravitational waves 
mentioned so far to the background gravitational radiation. We calculate f2 gw from Eqs. 
(13) and (14) for the highest luminosity case of L v = 6.5 x 10 52 erg s -1 . The resulting 
amplitude is presented in Figure 10. For the frequency greater than ~ 100 Hz, the signal 
is dominated by the contribution from matter motions, while at lower frequencies the signal 
is dominated by the one from the neutrinos (compare Figure 9). "Low rate" and "High 
rate"in the figure correspond to the cases of a = 0, 2 in Eq. (16), respectively, indicating 
the uncertainty of the core-collapse supernova rate as mentioned before. We terminated the 
simulation for the high luminosity case (see Table 1), seeing the shock wave propagate out 
of the iron core. For the frequencies of / <C (27rAt) _1 ~ 0.1 Hz(At/500 ms) -1 obtained 
from our simulations, we extended the spectra using the zero frequency limit (Turner 1978) 
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Time [ms] Time [ms] 

Fig. 7. — Effects of numerical resolutions on the waveforms. 60, 120 in the panel represents 
the number of the angular mesh points in the computational domain. Left and right panel 
shows the waveforms from the sum of the anisotropic neutrino emissions and only from the 
matter motions, respectively. After the saturation occurs (~ 200 ms), the difference of the 
resolution appears but not becomes larger than ~ 20% in the total amplitudes (left). 

as implemented in the study of Buonanno et al. (2005). "Upper bound" denotes the most 
optimistic estimation. In the case of L Ve = 6.5 x 10 52 erg s _1 , about 30 % of the available 
gravitational binding energy of the neutron star of 3 x 10 53 erg may be emitted during the 
simulation of 500 ms followed here. From Figure 2 showing the monotonic increase of the 
neutrino-originated amplitudes, we multiply the amplitude for the case of "High rate" by 
an enhancement factor of 3 and consider it as an upper bound. Then from Figure 10, we 
may say that the background radiation considered here could be larger for the frequency 
above ~ 1 Hz than the upper limit of the background radiation generated in the inflationary 
epoch (the horizontal line in the figure). However, the frequency range is just outside the 
sensitivity of the proposed detectors, such as DECIGO (Seto et al. 2001). 



4. Summary and Discussion 

We presented the results of numerical experiments, in which we studied how the as- 
phericities induced by the growth of the standing accretion shock instability could produce 
the gravitational waveforms in the postbounce phase of core-collapse supernovae. To ob- 
tain the neutrino-driven explosions, we parametrized the neutrino fluxes emitted from the 
central protoneutron star and approximated the neutrino transfer by the light-bulb scheme. 
By doing the spectrum analysis of the waveforms, we investigated the detectability of the 
signals from a single core-collapse supernova and the cosmological ones by the ground-based 
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Frequency [Hzl Frequency [Hzl 

Fig. 8. — Spectral distributions of the gravitational waves from the sum of anisotropic 
neutrino emissions and matter motions (thick line), and only from the matter motions (dotted 
line) for the case of L Ve = 6.5 (left) and 5.5 x 10 52 erg s -1 (right), respectively. It can be 
seen the dominance of the neutrino-originated gravitational waves at the lower frequency of 
~ 100 Hz. 

and space-based laser interferometers, respectively. Our main results can be summarized as 
follows. 

1. The amplitudes of the gravitational waves from the anisotropic neutrino emissions 
are larger up to two orders of magnitudes than the ones from the matter motions during 
the SASI operations. It is found that the wave amplitudes from the neutrinos show the 
monotonic increase with time, regardless of the neutrino luminosities from the protoneutron 
star. We point out that this feature can be understood by the specific nature of SASI, 
which makes the deformation of the shock waves of I = 1, 2 modes dominant, leading to the 
enhanced neutrino emissions in the regions close to the symmetry axis. In fact, we show 
that the amplitudes become larger when the growth of the SASI enters the nonlinear phase, 
in which the deformation of the shocks and the neutrino anisotropy become large. 

2. From the spectrum analysis of the waveforms, we find that the amplitudes from 
the anisotropic neutrino emissions are dominant over the ones from the matter motions at 
frequency < 100 Hz. The detection of such signals from a galactic supernova may be marginal 
for the currently running detector of the first LIGO and promising for the detectors in the 
next generation such as LCGT and the advanced LIGO. 

3. As for the background radiation, we indicate that the contribution of the gravitational 
signals considered here could be larger at frequency > 1 Hz than the primordial gravitational 
wave backgrounds generated in the inflationary epoch. Unfortunately, however, it is found 
that this frequency range is just outside of the sensitivity of the proposed detectors, such as 
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Fig. 9.— Detection limits of TAMA (Ando & the TAMA collaboration 2002), first LIGO 
(Thorne 1995), advanced LIGO (Weinstein 2002), and Large-scale Cryogenic Gravitational 
wave Telescope (LCGT) (LCGT Collaboration 1999) and the expected gravitational wave 
spectra. The source is assumed to be located at the distance of 10 kpc. 



DECIGO. 

We give a brief comparison with recently published models. The monotonic increase 
with time in the wave amplitudes of the neutrino-originated gravitational waves is consistent 
with the model sl5r of Miiller et al. (2004), in which the operation of the SASI was seen. 
However, the amplitudes here are typically larger (up to 1 order of magnitude). This should 
be mainly because the neutrino luminosity here is taken to be higher than the one obtained 
in Miiller et al. (2004), which is about ~ 2 x 10 52 erg/s during the SASI operation. As 
for the total gravitational-wave energy emission, typical values of the computed models 
here (~ 10~ lo M & c 2 , see Table 1) are one order magnitude smaller than the one in Miiller 
et al. (2004). This should be owing to the excision of the protoneutron star, by which 
the contribution from the high frequency domain of the energy emissions are eliminated in 
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this study. It should be noted that the larger oscillations of the protoneutron star in the 
postbounce phase (Burrows et al. 2006) and the resulting efficient gravitational emissions 
(Ott et al. 2006) cannot be treated in principle here. 

The simulation highlighted here is nothing but an idealized study for the physical un- 
derstanding of relation between the asphericities induced by the SASI and the resulting 
gravitational waves. Remembering the caveats about the assumptions of the artificially 
constructed initial condition, the fixed accretion rate, the absorbing boundary condition, 
and the fixed neutrino luminosity and energies, it is by no means definitive at all. Espe- 
cially, much better neutrino transfer is indispensable for more reliable calculations of the 
neutrino-originated gravitational wave, which we only considered the radial transport. One 
more major deficit is the axial symmetry assumed in the present two-dimensional (2D) sim- 
ulations. In three-dimensional environments, the pronounced dominance of / = 1,2 along 
the symmetry axis, which is a coordinate singularity in the 2D computations, may become 
weaker, owing to the additional spatial degree of freedom in the azimuthal (0) direction. In 
the 3D case, we think that the qualitative features of the plus mode waveform computed 
in the 2D case here will be unchanged, but quantitatively, we expect that the amplitudes 
become smaller owing to the reduced anisotropy along the symmetry axis. Thus the ampli- 
tudes calculated in this study could be an upper bound, in which the maximal anisotropy of 
the shock waves and thus neutrino emissions outside the neutrino sphere could be achieved. 
Furthermore we think that it is interesting to investigate the properties of the cross mode 
gravitational waves, which are of genuine 3D origin and could be possibly produced from the 
transfer of the / = 1, 2 modes to some modes with nonvanishing m in Yi m . This study is a 
prelude to the forthcoming 3D simulations to clarify those aspects, which will be presented 
elsewhere. 
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Fig. 10. — Cosmological gravitational wave background for the highest luminosity case of 
L Ue = 6.5 x 10 52 erg s _1 , taking into account the possible enhancement due to the un- 
certainty of the cosmic star formation rate labeled by "Low rate" and "High rate", and the 
limited simulation time labeled by "Upper bound" (see text for details), with some indicated 
sensitivity curves for future missions of LISA, the LIGO in the third generation (Buonanno 
et al. 2002), and practical/ultimate DECIGO (Seto et al. 2001). The horizontal line indi- 
cates a maximum version of gravitational wave backgrounds produced during the slow-roll 
inflation assuming T/S = 0.3 for the ratio of tensorial and scalar contributions to the cosmic 
microwave background radiation anisotropy (Turner 1997), concreted in Buonanno et al. 
(2005). 
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